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Abstract. This review demonstrates that neutrino emission from dense hadronic 
component in neutron stars is subject of strong modifications due to collective 
effects in the nuclear matter. With the most important in-medium processes incor- 
porated in the cooling code an overall agreement with available soft X ray data can 
be easily achieved. With these findings so called "standard" and "non-standard" 
cooling scenarios are replaced by one general "nuclear medium cooling scenario" 
which relates slow and rapid neutron star coolings to the star masses (interior den- 
sities). In-medium effects take important part also at early hot stage of neutron 
star evolution decreasing the neutrino opacity for less massive and increasing for 
more massive neutron stars. A formalism for calculation of neutrino radiation from 
nuclear matter is presented that treats on equal footing one-nucleon and multiple- 
nucleon processes as well as reactions with resonance bosons and condensates. Cool- 
ing history of neutron stars with quark cores is also discussed. 

1 Introduction 

The EINSTEIN, EXOSAT and ROSAT observatories have measured surface 
temperatures of certain neutron stars (NS) and put upper limits on the sur- 
face temperatures of some other NS ( cf. [1-3] and further references therein). 
The data for some supernova remnants indicate rather slow cooling, while the 
data for several pulsars point to an essentially more rapid cooling. 

Physics of NS cooling is based on a number of ingredients, among which 
the neutrino emissivity of the high density hadronic matter in the star core 
plays a crucial role. Neutron star temperatures are such that, except first 
minutes-hours, neutrinos/antineutrinos radiate energy directly from the star 
without subsequent collisions, since Ai,,A p 3> -R, where A„,A P are the neu- 
trino and antineutrino mean free paths and R is star radius. In the so called 
"standard scenario" of the NS cooling (scenario for slow cooling) the most im- 
portant channel up to temperatures T ~ 10 9 K belongs to the modified Urea 
(MU) process nn — > npeD. First estimates of its emissivity were done in 
[4,5]. References [6,7] recalculated the emissivity of this process in the model, 
where the nucleon-nucleon (NN) interaction was treated with the help of 
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slightly modified free one-pion exchange (FOPE). This important result for 
the emissivity, e v [FOPE] , was proved to be by an order of magnitude larger 
than the previously obtained one. Namely the value e„[FOPE] was used in 
various computer simulations resulting in the "standard scenario" of the cool- 
ing, e.g. cf. [8-10]. Besides the MU process, in the framework of the "standard 
scenario " numerical codes included also processes of the nucleon (neutron and 
proton) bremsstrahlung (NB) nn^> nnvv and np — > npvv, which lead to a 
smaller contribution to the emissivity then the MU, cf. [11,6]. Medium effects 
enter the above two-nucleon (MU and NB) rates mainly through the effec- 
tive mass of the nucleons which has a smooth density dependence. Therefore 
within FOPE model the density dependence of the reaction rates is rather 
weak and the neutrino radiation from a NS depends only very weakly on its 
mass. This is the reason why the "standard scenario" based on the result [6], 
though complying well with several slowly cooling pulsars, fails to explain the 
data of the more rapidly cooling ones. Also "standard scenario " included pro- 
cesses contributing to the emissivity in the NS crust which become important 
at a lower temperature. 

The non-standard scenario included so called exotica, associated with 
different types of direct Urca-like processes, i.e. the pion Urea (PU) [12] 
and kaon Urea (KU) [13,14] /3-decay processes and direct Urea (DU) on 
nucleons and hyperons [15] possible only in sufficiently dense interiors of 
rather massive NS. The main difference in the cooling efficiency driven by the 
DU-like processes on one hand and the MU and NB processes on the other 
hand lies in the rather different phase spaces associated with these reactions. 
In the MU and NB case the available phase space is that of a two-fermion 
origin, while in the pion (kaon) /3-decay and DU on nucleons and hyperons 
it is that of a one-fermion origin. Critical density of pion condensation in NS 
matter is g C7r ~ (1 -j- 3) go depending on the type of condensation (neutral or 
charged) and the model, see [16-19]. Critical density of kaon condensation 
is q c k — (2 -j- 6) go depending on the type (K ~ or K°, S or P wave) and 
the model, see [20,21]. Critical density for the DU process is g c u ~ (2 6) go 
depending on the model for the equation of state (EQS), see [15,18]. Recent 
calculations [19] estimated critical density of neutral pion condensation as 
2.5go and for the charged one as 1.7go> whereas variational calculations [18] 
argued even for smaller critical densities (~ 1.3 go for 7r° condensation). On 
the other hand the EQS of [18] allows for DU process only at g > 5go- 

There is no bridge between "standard" and "non-standard" scenarios due 
to complete ignorance of in-medium modifications of NN interaction which 
allows for strong polarization of the soft modes (like virtual dressed pion 
and kaon modes serving a part of in- medium baryon-baryon interaction). 
Only due to enhancement of a medium polarization with the baryon density 
pion condensate may appear and it seems thereby quite inconsistent to ignore 
the softening effects for g < g C7T , and suddenly switch on the condensate for 

g > gen- 
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Now let us basing on the results [22-25,16,26] briefly discuss a general 
"nuclear medium cooling scenario " which treats obvious caveats of two men- 
tioned above scenarios. First of all one observes [22] that in the nuclear mat- 
ter many new reaction channels are opened up compared to the vacuum 
processes. Standard Feynman technique fails to calculate in-medium reac- 
tion rates if the particle widths are important since there are no free particle 
asymptotic states in matter. Then summation of all perturbative Feynman 
diagrams where free Green functions are replaced by the in-medium ones 
leads to a double counting due to multiple repetitions of some processes (for 
an extensive discussion of this defect see [27]). This calls a formalism dealing 
with closed diagrams (integrated over all possible in-medium particle states) 
with full non-equilibrium Green functions. Such a formalism was elaborated 
in [23,24] first within quasiparticle approximation (QPA) for nucleons and was 
called in [24] "optic theorem formalism (OTF) in non- equilibrium diagram 
technique". It was demonstrated that standard calculation of the rates via 
squared reaction matrix elements and calculation using OTF coincide within 
QPA picture for the fermions. In [28] the formalism was generalized to in- 
clude arbitrary particle widths effects. The latter formalism treats on equal 
footing one-nucleon and multiple-nucleon processes as well as resonance re- 
action contributions of the boson origin, as processes with participation of 
zero sounds and reactions on the boson condensates. Each diagram in the 
series with full Green functions is free from the infrared divergences. Both, 
the correct quasiparticle (QP) and quasiclassical limits are recovered. 

Except for very early stage of NS evolution (minutes - hours) typical 
averaged lepton energy (> T) is larger then the nucleon particle width TV ~ 
T 2 /efn and the nucleons can be treated within the QPA. This observation 
much simplifies consideration since one can use an intuitive way of separation 
of the processes according to the available phase space. The one-nucleon 
processes have the largest emissivity (if they are not forbidden by energy- 
momentum conservations), then two-nucleon processes come into play, etc. 

In the temperature interval T c < T < T opac (T c is typical tempera- 
ture for the nucleon pairing and T opac is typical temperature at which neu- 
trino/antineutrino mean free path A^/Ap is approximately equal to the star 
radius R) the neutrino emission is dominated by the medium modified Urea 
(MMU) and medium nucleon bremsstrahlung (MNB) processes if one-nucleon 
reactions like DU, PU and KU are forbidden, as it is the case for g < 
QcU i Qcix, QcK- Corresponding diagrams for MMU process are schematically 
shown in Fig. 1. References [22-25,16] considered NN interaction within 
Fermi liquid Landau-Migdal approach. They incorporated the softening of 
the medium one-pion exchange (MOPE) mode, other medium polarization 
effects, like nucleon-nucleon correlations in the vertices, renormalization of 
the local part of NN interaction by the loops, as well as the possibility of the 
neutrino emission from the intermediate reaction states and resonance DU- 
like reactions going on zero sounds and the boson condensates. References 
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Fig. 1. Antincutrino emission from a nucleon leg (left graph) and from intermediate 
scattering states (right) in MMU process. Full dot includes weak coupling vertex 
rcnormalization. 



[23,25,16] have demonstrated that for g > go second diagram of Fig. 1 gives 
the main contribution to the emissivity of MMU process rather than the first 
one which contribution has been earlier evaluated in the framework of FOPE 
model in [6]. This fact essentially modifies the absolute value as well as the 
density dependence of the nn — > npeP process rate which becomes to be very 
sharp. Thereby, for stars of masses larger than the solar mass the resulting 
emissivities were proved to be substantially larger than those values calcu- 
lated in FOPE model. With increase of the star mass (central density) pion 
mode continues to soften and MMU and MNB rates still increase. At g > g C7r 
pion condensation begins to contribute. Actually, the condensate droplets 
may be exist already at a smaller density if the mixed phase is appeared, 
as suggested in [29]. At T > T me i t , where T me i t is the melting temperature, 
roughly <~ several MeV, the mixed phase is in liquid state and PU processes 
on independent condensate droplets are possible. At T < T me i t condensate 
droplets are placed in a crystalline lattice that substantially suppresses cor- 
responding neutrino processes. 

Reference [12] considered the reaction channel n —> pir~ eP, whereas [22,23] 
included other possible pion 7r + , n ± , ir° condensate processes with charged 
and neutral currents (e.g., like mr® — ► peP, mr+ — > pvv and mr® — > nvP) 
as well as resonance reactions going on zero sounds which are also possible 
at g < g clx . Due to NN correlations all pion condensate rates are signifi- 
cantly suppressed (by factors ~ 10 - 100 compared to first estimate [12], see 
[30,22,23,31]. At g ~ g c1T both MMU and PU processes are of the same order 
of magnitude [23] demonstrating a smooth transition to higher densities (star 
masses) being absent in the "standard" and "non-standard" scenarios. 

For T < T c the reactions of neutrino pair radiation from supcrfluid 
nucleon pair breaking and formation (NPBF) shown in Fig. 2 become to 
be dominant processes. The neutron pair breaking and formation (nPBF) 
process for the case of the ISq pairing was first calculated in [32] using 
standard Bogolyubov technique. Later this process was independently cal- 
culated in [24,25] as demonstration of efficiency of OTF within the closed 
non-equilibrium diagram technique developed there. Moreover [24] calculated 
emissivity of the corresponding process on proton (pPFB) taking into account 
strong coupling pvv vertex renormalization (see first diagram (19) below). It 
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Fig. 2. Neutrino-antineutrino emission from Cooper pair-breaking (left graph) and 
pair- formation processes (right graph). 

results in one-two order of magnitude enhancement of pPFB cmissivity (for 
q <~ (1.5 -j- 3)po) compared to that would be estimated with the vacuum ver- 
tex, leading to that both nPFB and pPFB emissivities can be of the same 
order of magnitude. Emissivities of NPBF processes have the same suppres- 
sion factor - exp(-2Z\/T) as MU, NB, MMU and MNB at T < T c but 
compared to the latter the NPBF processes have a large one nucleon phase 
space volume. Ref. [24] and then [16] also sketched how to incorporate 3P2 
pairing (the developed formalism easily allows to do it) indicating that expo- 
nential suppression of the specific heat and the emissivity is to be replaced 
by only a power law suppression when the gap vanishes at the Fermi sphere 
poles. This is precisely the case for the 3P2G mj \— 2) pairing, where mj is 
projection of the total pair momentum onto quantization axis. This idea was 
then worked out in [33-35] , without reference. 

I would also like to do a historic remark relating to estimation of the 
NPBF processes since in some works (e.g., see [36,35]) was expressed a sur- 
prise why these processes were not on a market during many years. First 
work [32], although found correct analytic expression for 1S pairing of neu- 
trons, numerically underestimated the emissivity by an order of magnitude. 
Also there was no statement on the dominance of the process in the cool- 
ing scenario (i.e. over the MU). Asymptotic behaviour of the emissivity 
e[nPBF] - 10 20 T 9 7 exp(-2Z\/T) for T < A, T 9 = T/10 9 K, as follows from 
expression (lb) of [32] and from their rough asymptotic estimate of the inte- 
gral (see below (13b)), shows up nor a large one nucleon phase space factor 
(<~ 10 28 ) nor appropriate temperature behaviour. Namely this underestima- 
tion of the rate in [32] (we again point out that analytic expression (la) 
is correct) and absence of mentioning of possible dominance of the process 
over MU, became the reason that this important result was overlooked dur- 
ing many years. Reference [24] overlooked sign of anomalous diagram for 
S pairing and included oc g\ term (g^ is axial-vector coupling constant) 
which should absent for S pairing giving, nevertheless, the main contribution 
for the 3P2 pairing case. A reasonable numerical estimate of the emissiv- 
ity was presented valid for both S and P pairings (except | mj | = 2 case) 
including NN correlation effects into consideration. Uncertainty of this esti- 
mate is given by a factor (0.5 -j- 2) which is allowed by variation of not too 
well known correlation factors. Namely this estimate was then incorporated 
within the cooling code in [26]. Correct asymptotic behaviour of the emis- 
sivity is e[nPFB, pPFB] - 10 28 (Z\/McV) 7 (T/Z\) 1 / 2 cxp(-2Z\/T) for T < A, 



6 D.N. Voskresensky 



that showed up a huge one-nucleon phase space factor and very moderate T 
dependence of the pre-factor. Thereby the value of the rate was related in 
[24] to the value of the pairing gap. The possibility of the dominant role of 
this process (even compared to enhanced MMU and PU rather than only to 
MU) was unambiguously stressed. Unfortunately [24] had a number of obvi- 
ous misprints which were partially corrected in subsequent papers. Although 
it does not excuse the authors these misprints can be easily treated by an 
attentive reader. Importance of NPFB processes was then once more stressed 
in review [16]. Reference [26] was the first that quoted the previous result 
[32]. It incorporated most important in- medium effects in the cooling code, 
among them nPBF and pPBF as equally important processes. The latter 
process was then rediscovered in [34] as negligible one due to ignorance there 
of correlations in weak interaction vertices. Work [36] supported conclusion 
of [26] on importance of NPBF processes as governing the cooling scenario 
at T < T c . 

The medium modifications of all the above mentioned rates result in a 
pronounced density dependence (for NPBF processes mainly via dependence 
of the pairing gaps on the density and dependence on the NN correlation 
factors), which links the cooling behavior of a neutron star decisively to its 
mass [22,23,16,26]. As the result, the above mentioned medium modifications 
lead to a more rapid cooling than obtained in the "standard scenario ". Hence 
they provide a possible explanation for the observed deviations of some of the 
pulsar temperatures from the "standard" cooling. Particularly, they provide 
a smooth transition from "standard" to "non-standard" cooling for increasing 
central star densities, i.e., star masses. Thus by means of taking into account 
of most important in-medium effects in the reaction rates one is indeed able 
to achieve an appropriate agreement with both the high as well as the low 
observed pulsar temperatures that leads to the new "nuclear medium cooling 
scenario". Using a collection of modern EQS for nuclear matter, which cov- 
ered both relativistic as well as non-relativistic models, [26] also has demon- 
strated a relative robustness of these in-medium cooling mechanisms against 
variations in the EQS of dense NS matter (for EQS that allow for a wide 
dense hadronic region in NS). 

At initial stage (T > T opac ) the newly formed hot NS is opaque for neu- 
trinos/antineutrinos. Within FOPE model the value T opac was estimated in 
[6]. Elastic scatterings were included in [37,38] and pion condensation effect 
on the opacity was discussed in [39] . Medium effects dramatically affect the 
neutrino/antineutrino mean free paths, since A„(p) oc 1/ | M | 2 , where M 
is the reaction matrix element. Thereby, one-nucleon elastic scattering pro- 
cesses, like Nv — > Nv, for energy and momentum transfer lu < qvFN are 
suppressed by NN correlations [24,16,40,41]. Neutrino/antineutrino absorp- 
tion in two-nucleon MMU and MNB processes is substantially increased with 
the density (since | M | 2 for MMU and MNB processes increases with the 
density) [23,24,16]. Thus more massive NS are opaque for neutrinos up to 
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lower temperatures that also results in a delay of neutrino pulse. Within the 
QPA for the nucleons the value T opac was estimated with taking into ac- 
count of medium effects in [23,16]. References [42,43,16] considered possible 
consequences of such a delay for supernova explosions. On the other hand, 
at T > (1-7- 2)MeV one should take care of the neutrino/antineutrino ra- 
diation in multiple NN scatterings (Landau-Pomeranchuk-Migdal (LPM) 
effect) when averaged neutrino-antineutrino energy, uj„ p ~ several T, be- 
comes to be smaller than the nucleon width [28]. Numerical evaluations 
of Jjv in application to MNB processes were done in [44] and [45] within the 
Bruckner scheme and the Bethe-Salpeter equation, respectively. The LPM 
effect suppresses the rates of the neutrino elastic scattering processes on nu- 
cleons and also it suppresses MNB rates. For NS of rather low mass ( Mq ) 
the suppression of the rates of neutral current processes due to the multiple 
collision coherence effect prevails over the enhancement due to the pion soft- 
ening, and for sufficiently massive NS (> 1.4Mq) the enhancement prevails 
the suppression. MMU cmissivity remains to be almost unaffected by the 
LPM effect since averaged De energy ~ pp e is rather large (3> -Tat). 

Besides pion and kaon condensates, under current discussion is the pos- 
sibility of the presence of the quark matter at sufficiently high baryon den- 
sities. Thereby, these phases may exist in most massive NS. Also a partic- 
ular discussion is devoted to the possibility of the existence of less mas- 
sive, stable or metastable, dense nuclei-stars glued by the condensates or by 
strange quarks, cf. [46,47,16,48,49]. Uncertainty in these predictions relate 
to the poorly known EQS at high baryon densities. Several possible types of 
stars differing in their values of the bag constant B were discussed: ordinary 
NS without any quark core, hybrid neutron stars (HNS) with quark matter 
present only in their deep interiors (for somewhat intermediate values of B), 
NS with large quark core (QCNS), a narrow hadron layer and a crust typical 
for a NS, and quark stars (QS) with a tiny crust of normal matter and with 
no crust (both for low B). If these objects are produced in supernova explo- 
sions within ordinary mechanism of blowing off the mantle, they are rather 
massive HNS. If an extra support for the blowing off the matter arises, there 
may appear less massive objects like QCNS, QS and may be even objects of 
arbitrary size. 

Especial interest for the discussion of cooling of HNS, QCNS, and QS 
is motivated by recent works [50,51] which demonstrated possibility of di- 
quark condensates characterized by large pairing gaps (A q 100 MeV) in 
dense quark matter. The two-flavor (2SC) and the three-flavor (3SC) color 
superconducting phases allow for unpaired quarks of one color whereas in 
the color-flavor locking (CFL) phase all the quarks are paired. Presence of 
large quark gap may significantly affect the cooling history of HS, QCNS and 
QS. Thus it is interesting to confront this possibility to the observational 
data. Cooling of QCNS, and QS was discussed in [52,53] and the case of HNS 
was considered in [54,53]. In difference with [52,54] more recent work [53] 
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incorporated the heat transport being important during a longer period of 
time compared to that can be in absence of the color superconducting quark 
matter. 

The paper is organized as follows. Sect. 2 discusses basic ideas of the 
Fermi liquid approach to description of nuclear matter. The NN interaction 
amplitude is constructed with an explicit treatment of long-ranged soft pion 
mode and vertex rcnormalizations due to NN correlations. The meaning of 
the pion softening effect is clarified and a comparison of MOPE and FOPE 
models is given. Also renormalization of the weak interaction in NS matter is 
performed. Sect. 3 discusses the cooling of NS at T < T opac . Comparison of 
cmissivitics of MMU and MU processes shows a significant enhancement of in- 
medium rates. Then we discuss DU-like processes and demonstrate medium 
effect due to vertex rcnormalizations. The role of in-medium mechanisms in 
the cooling evolution of NS is then demonstrated within a realistic cooling 
code. Next we consider influence of in-medium effects on the neutrino mean 
free path at initial stage of NS cooling. Essential role of multiple NN colli- 
sions is discussed. Sect. 4 presents OTF in non-equilibrium closed diagram 
technique in the framework of QPA for the nucleons and also beyond the QPA 
incorporating genuine particle width effects. In the last Sect. 5 following [53] 
we review recent results on the cooling of HNS. 

2 Nuclear Fermi liquid description 

2.1 NN interaction. Separation of hard and soft modes 

At temperatures of our interest (T <C £f„) neutrons are only slightly excited 
above their Fermi sea and all the processes occur in a narrow vicinity of 
Epn- In such a situation Fermi liquid approach seems to be the most efficient 
one. Within this approach the long-ranged diagrams are treated explicitly 
whereas short-scale diagrams are supposed to be the local quantities given 
by phenomenological so called Landau-Migdal (LM) parameters. Thus using 
argumentation of Fermi liquid theory [55-57,16] the retarded NN interaction 
amplitude is presented as follows (see also [58]) 




(1) 



where 





The solid line presents the nucleon, whereas double-line, the A isobar. The 
double- wavy line corresponds to the exchange of the free pion with inclusion 
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of the contributions of the residual S wave ttNN interaction and tttt scat- 
tering, i.e. the residual irreducible interaction to the nucleon particle-holes 
and delta-nucleon holes. The full particle-hole, delta-nucleon hole and pion 
irreducible block (first block in (2)) is by its construction essentially more 
local then contributions given by explicitly presented graphs. Thereby, it is 
parameterized with the help of the LM parameters. In the standard Landau 
Fermi liquid theory fermions are supposed to be at their Fermi surface and 
the Landau parameters are further expanded in Lcgcndre polynomials in the 
angle between fermionic momenta. Luckily only zero and first harmonics en- 
ter physical quantities. The momentum dependence of the residual part of 
nuclear forces is expected to be not as pronounced and one can avoid this 
expansion. Then these parameters, i.e. f nn , f np and g nn , g np in scalar and 
spin channels respectively are considered as slightly momentum dependent 
quantities. In principle, they should be calculated as functions of the density, 
neutron and proton concentrations, energy and momentum but, simplifying, 
one can extract them from analysis of experimental data. 

The part of interaction involving A isobar is analogously constructed 




The main part of the N A interaction is due to the pion exchange. Although 
information on local part of the N A interaction is rather scarce, one can 
conclude [16,19] that the corresponding LM parameters are essentially smaller 
then those for NN interaction. Besides, at small transferred energies uj <C m x 
the Z\-nucleon hole contribution is a smooth function of uj and k in difference 
with the nucleon-nucleon hole (NN -1 ) contribution. Therefore and also for 
simplicity we neglect the first graph in r.h.s. of (3). 

Straightforward resummation of (1) in neutral channel yields [24,16] 

P 

r*p = = Co n + 3>i ■ v 2 ) + H n T*(ct x ■ fc)(<x 2 • fc), (4) 

a (3 

T% = faplifafl), Z* n = g n „7(g„„), = g npl (g nn ), a, (3 = (n,p), 

Tun = J 2 (gnn)D% h T« - -7pp7(ffn„)^o , T* = 7p 2 p £«o, 

1-\x) = l-2xC A* n , 7 PP = (1 - 4 ff Gv4«J 7 ( 5 „„), (5) 

fnn = fpp = / + /', fnp = f - f, g nn = g PP = .9 + and 9n P = 9 - 9' , di- 
mensional normalization factor is usually taken to be Co = tt 2 /[itinPf(qo)] — 
300 MeV- fm 3 ~ 0.77m" 2 , D^ a is the full retarded Green function of tt°, A afj 
is the corresponding NN -1 loop (without spin degeneracy factor 2) 
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A ^ = O . ^n(w q) * (47T 2 )" 1 - 2vf„) , (6) 

a- 1 

j4„ rl ~ — m^pp n (2n 2 )~ 1 , for ui <C qv Fni q <C 2p F „, and we for simplicity 

neglect proton hole contributions due to a small concentration of protons. 

Resummation of (1) in the charged channel yields 
n n 

r* = ^ = Co (X« + Z% <7l • <X 2 ) + /^T* (<Ti • fc) (<T 2 • fc) , (7) 

p p 

T% = 2/'7(/'), Z r f p = 2^7(5'). = tVP* , (8) 
7~ 1 W = l-4xC A% p . 

The LM parameters are rather unknown for isospin asymmetric nuclear mat- 
ter and for g > g . Although some evaluations of these quantities have been 
done, much work is still needed to get convincing results. Therefore for es- 
timates we will use the values extracted from atomic nucleus experiments. 
Using argumentation of a relative locality of these quantities we will suppose 
the LM parameters to be independent on the density for g > qq. One then 
can expect that the most uncertain will be the value of the scalar constant 
/ due to essential role of the medium-heavy a meson in this channel. But 
this parameter does not enter the tensor force channel being most important 
in our case. Unfortunately there are also essential uncertainties in numerical 
values of some of the LM parameters even for atomic nuclei. These uncertain- 
ties are, mainly, due to attempts to get the best fit of experimental data in 
each concrete case slightly modifying parameterization used for the residual 
part of the NN interaction. E.g., calculations [56] gave / ~ 0.25, /' ~ 1, 
g ~ 0.5, g' ~ 1 whereas [59-61], including QP renormalization pre-factors, 
derived the values / ~ 0, /' ~ 0.5 + 0.6, g ~ 0.05 ± 0.1, g 1 ~ 1.1 ± 0.1. 

Typical energies and momenta entering NN interaction of our interest 
are w ~ and k ~ PFn- Then rough estimation yields 7(<7rm,w — 0, k ~ 
PFm Q — Qo) — 0.35 -j- 0.45. For w = k ~ T typical for the weak processes 
with participation of vv one has -f^ 1 (g nn , w ~ k ~ T, g = g ) ~ 0.8 -j- 0.9. 



2.2 Virtual Pion Mode 

Resummation of diagrams yields the following Dyson equation for pions 



+ 



77£ 



(9) 
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The irNA full-dot-vertex includes a background correction due to the pres- 
ence of the higher resonances, Ilf} es is the residual retarded pion self-energy 
that includes the contribution of all the diagrams which are not presented ex- 
plicitly in (9), as S wave ir TV TV and tttt scatterings (included by double- wavy 
line in (2)). The full vertex takes into account TV TV correlations 




(10) 



Due to that the nucleon particle- hole part of II ^ is oc j(g nn ) an d the nucleon 
particle-hole part of II 7r ± is oc "f(g') ■ The value of the TV TV interaction in the 
pion channel is determined by the full pion propagator at small u) and k ~ 
PFn, i-e. by the quantity uj 2 (k) — —(D^)^ 1 ^ = 0, k,^ % ). Typical momenta 
of our interest are k ~ PFn- Indeed the momenta entering the TV TV interaction 
in MU and MMU processes are k = pm, the momenta governing the MNB 
are k = ko [23] where the value k = fc ~ (0.9 -j- l)pFn corresponds to the 
minimum of uj 2 (k). The quantity ui = u>(ko) has the meaning of the effective 
pion gap. It is different for tt° and for ir^ since neutral and charged channels 
are characterized by different diagrams permitted by charge conservation, 
thus also depending on the value of the pion chemical potential, =/= fi^- 
0, = 0. For T <C £i?„,£i? p , one has /i^- = fj, e = ep n — Ef p , as follows 
from equilibrium conditions for the reactions n — > pn~ and n — > pev. 

Change of the sign of uj 2 symbolizes the pion condensate phase transition. 
Typical density behaviour of u 2 is shown in Fig. 3. At g < (0.5 -j- 0.7)qq, one 




Fig. 3. Effective pion gap (for fi^o = 0) versus baryon density, see [16]. 



has uj 2 = m 2 , — /j, 2 . For such densities the value u> 2 (pFn) essentially deviates 
from m 2 — fi 2 tending to m 2 +p% n — ^ in small density limit. 

At the critical point of the pion condensation (g = g C7r ) the value uj 2 with 
artificially switched off tttt fluctuations (dashed line in Fig. 3) changes its sign. 
In reality tttt fluctuations are significant in the vicinity of the critical point 
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and there occurs the first-order phase transition to the inhomogeneous pion- 
condensate state [62-64]. Thereby there are two branches (solid curves in Fig. 
3) with positive and respectively negative values for uj 2 . Calculations of [64] 
demonstrated that at g > g C7T the free energy of the state with uj 2 > 0, where 
the pion mean field is zero, becomes larger than that of the corresponding 
state with uj 2 < and a finite mean field. Therefore at g > g C7r the state 
with uj 2 > is metastable and the state with uj 2 < and the pion mean field 
(fn ^ becomes the ground state. 

The quantity uj 2 demonstrates how much the virtual (particle-hole) mode 
with pion quantum numbers is softened at given density. The ratio a = 
Z) 7r [med.]/-D 7r [vac.] ~ 6 for g = go, uj = 0, k = pfn and for isospin symmetric 
nuclear matter. However this essential so called "pion softening" [57] does not 
significantly enhance the NN scattering cross section due to a simultaneous 
essential suppression of the ttNN vertex by NN correlations. Indeed, the 
ratio of the NN cross sections calculated with FOPE and MOPE is 



and for g = g we have R & 1, whereas for g = 2g we already get R ~ 10. 

As follows from numerical estimates of different 7 factors entering (4) and 
(7), the main contribution to NN interaction for g > go is given by MOPE 



whether this channel (T oc (eri ■ fc)(<r 2 ■ k)) of the reaction is not forbidden or 
suppressed by some specific reasons like symmetry, small momentum transfer, 
etc. 

The g meson contribution to NN interaction is partially included in g a p, 
other part contributing to T and T is minor (oc uj 2 /m e 2 ). Indeed, using that 
(eri x fc)(er 2 x k) = fc 2 cr 1 cr2 — (cr 1 k) ■ (er 2 fc) the p exchange can be cast as 



for p <~ m\ under consideration; 7 factor is the same as for the corresponding 
pion (charged or neutral). First term is supposed to be included in phe- 
nomenological value of the corresponding Landau-Migdal parameter leading 
to its momentum dependence. Due to that the value g' gets a 30% decrease 
(rather then an increase discussed in some works) with the momentum at 
k = pfn compared to the corresponding value at k = 0, see [65]. Thus 
one can think that g'(0) > g'(j>FN)- Second term can be dropped as small 
( £ 1/30 of MOPE contribution at p = p ). 



R = 



q[MOPE] 7 4 (g>-0,fc-p FW )K+p 2 Fjv ) 2 
cr[FOPE] to 4 (pfn) 



(11) 



I 




(12) 
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Another worry [66,67] was expressed in connection with quasielastic po- 
larization transfer experiment at LAMPF, EMC experiment and the Drell- 
Yan experiment at Fermi Lab which did not observe an expected in several 
works a pronounced pion excess in nuclei. First optimistic estimates demon- 
strated the pion excess to be at (15-30)% level that was ruled out by different 
models which analyzed mentioned experimental results including artificially 
enhanced pion contribution and suppressing other important contributions. 
However it may well be that only in a very narrow vicinity of g C7r (far beyond 
go) pion fluctuations grow at T — and only at T ^ pion excess must 
essentially grow in substantially wider density interval [63,64]. 

Let us show that some effect should be although it may be numerically 
rather small. Let estimate a contribution of the virtual pion sea (region of 
small energies u) < kv F ) to the sum rule. It is given by 



f kVF dw 2 

S(kv F ) = / 2wAr — = — 
Jo 2lT 



kvF — 



2 (fc) 



arctan 



( pkv F \ 



« 1,(14) 



is the pion spectral function, 5(oo) = 1, f3 



'■N 



is the non- 



relativistic effective mass given by the relation (de(p)/d p)\ PF = pf/iti^. The 
quantity S(kvp) having no singularity even at u> 2 (ko) is, thus, rather 
insensitive to the pion softening effect. Therefore mentioned experimental 
analysis might be not as critical to the value w 2 (fco). 

The physical reason of a contribution of virtual pions at w < kvp is 
clear. It is the well known Landau damping associated with possibility of the 
virtual pion decay to the nucleon particle and hole, i.e. just with the Pauli 
blocking. As known from irN scattering experiments, pions interact with 
nucleons. Thus must be a contribution of virtual pions to the sum-rule for 
pion spectral function. There is no other way out. Therefore on a question 
"Where virtual pions are?" we would still suggest an old fashion answer: 
" They are hidden inside the matter" . The question which remains is only 
"What is the actual value of the pion excess?" To properly answer it one needs 
besides the pionic contribution to elaborate all relevant specific contributions 
related to the phenomenon analyzed in the given concrete experiment. E.g., 
quasielastic polarization transfer data are reasonably fitted with a slightly 
increased value of the Landau-Migdal parameter g' and with inclusion of S 
wave repulsion at small energies [16,68,69], EMC experiment is reproduced 
with inclusion of the soft pion mode and other relevant effects [70] . 

Thus instead of FOPE+g exchange, as the model of NN interaction which 
has been used in [6] in calculation of the emissivities of the two-nucleon 
reactions within the "standard scenario" of NS cooling, one should use the 
full NN interaction given by (4) and (7) or, simplifying, approximated by its 
MOPE part. 



14 D.N. Voskresensky 



2.3 Renormalization of the weak interaction. 

The full weak coupling vertex that takes into account NN correlations is 
determined by (10) where now the wavy line should be replaced by the lcpton 
pair. Thus for the vertex of our interest, N\ — > N^lv, we obtain [24,16] 

V fi = ^=m')lo-9AW)l<r], (15) 

for the /3 decay and 

G G 

Vnn = --^J^ hUnn)h ~ gAl{gnn)lv] , V^, = -^-j= [Hp p l - gAlpplv] (16) 

n PP = c v - 2f np j(f nn )C A nn , jp P = (1 - 4gC A rm ) j(g nn ), (17) 

for processes on the neutral currents N1N2 — > N1N2VV, V pp = V p ^ + V^, 
G ~ 1.17 • 10~ 5 GcV~ 2 is the Fermi weak coupling constant, cy = 1 — 
4sin 2 #i4'> sin 2 9w — 0.23, 5,4 ~ 1.26 is the axial-vector coupling constant, 
and lp = u{qi)^p{l — 75)11(92) is the lepton current. The pion contribution 
~ q 2 is small for typical | q |~ T or pFe, and for simplicity is omitted. In 
medium the value of gA (i.e. 3^) slightly decreases with the density that can 
be easily incorporated, cf. the Brown-Rho scaling idea [71]. Please notice 
that with a decrease of g* Al in particular, the value g cw increases remaining 
however to be finite (~ 2qq according to [72] where this decrease of g* A up to 1 
was discussed) due to attractive contribution of the A isobar in (9), whereas 
one would expect g C7I — > 00 for g* A — * 1 ignoring A contribution. 

The 7 factors renormalize the corresponding vacuum vertices. These fac- 
tors are essentially different for different processes involved. The matrix el- 
ements of the neutrino/ antineutrino scattering processes Nv — > Nv and of 
MNB behave differently in dependence on the energy-momentum transfer 

and whether N — n or N = p in the weak coupling vertex. Vertices 
v v v 



(18) 



N N N 

are modified by the correlation factors (5) and (8). For N — n these are 

j(gnn,w,q) and j(fnn,v, q) leading to an enhancement of the cross sections 

for lu > qvFn and to a suppression for lu < qvFn- Renormalization of the 

proton vertex (vector part of V p ^ + V^ p ) is governed by the processes [24,73] 
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being forbidden in vacuum. For the systems with ISo proton-proton pairing, 
oc g A contribution to the squared matrix element (see (16)) is compensated by 
the corresponding contribution of the diagram with anomalous Green func- 
tions of protons. The vector current term is oc Cy in vacuum whereas it is 
oc Kp p in medium (by first graph (19)). Thereby the corresponding vertices 
with participation of proton are enhanced in medium compared to small vac- 
uum value (oc c v ~ 0.006) leading to enhancement of the cross sections, up to 
<~ 104- 10 2 times for 1.543g»o depending on parameter choice. It does not con- 
tradict to the statement of ref. [24] that correlations are rather suppressed 
in the weak interaction vertices at g < g . Enhancement with the density 
comes from estimation (6) that directly follows from equation (2.30) of [24]. 
Also enhancement factor (up to ~ 10 2 ) comes from the virtual in-medium 
photon (j m ) whose propagator contains 1/m 2 <~ 1/e 2 , where m 7 is the ef- 
fective spectrum gap, that compensates small e 2 factor from electromagnetic 
vertices, see [73]. It is included by replacement Vj£ — > V pp . Other processes 
permitted in intermediate states like processes with pp' 1 and with the pion 
are suppressed by a small proton density and by q 2 ~ T 2 pre-factors, respec- 
tively. First diagram (19) was considered in [24], where the pPBF process was 
suggested, and then in [16,26], and it was shown that nPBF and pPBF pro- 
cesses may give contributions of the same order of magnitude. Finally, with 
electron and nucleon correlations included we recover this statement and nu- 
merical estimate [24,26] (in spite of mentioned misprints in the result [24]). 
Several subsequent papers [34,74,35] rediscovered pPBF process ignoring nu- 
cleon and electron correlation effects. (Although the authors were personally 
informed, they continue to insist [75] on incorrect treatment, giving prior- 
ity to their result.) Contribution of second diagram for pPFB process was 
recently incorporated in [76]. 

Paper [27] gives another example demonstrating that, although the vac- 
uum branching ratio of the kaon decays is r(K~ — > e~ + v e )/r(K~ — ► 
p~ + v^) w 2.5 x 10 -5 , in medium (due to Ap^ 1 decays of virtual K~) it 
becomes to be of order of unit. Thus we again sec that in dependence of what 
reaction channel is considered in-medium effects may or strongly enhance the 
reaction rates or substantially suppress them. Ignorance of these effects may 
lead to misleading results where one struggling for numerical factors ~ 1 may 
loose by orders of magnitude larger factors. 

2.4 Inconsistencies of FOPE model 

Since FOPE model became the base of the "standard scenario" for cooling 
simulations we would like first to demonstrate principal inconsistencies of the 
model for the description of interactions in dense (g Z go) baryon medium. 
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The only diagram in FOPE model which contributes to the MU and NB is 

(20) 




o 



firN 

Dots symbolize FOPE. This is first available Born approximation diagram, 
i.e. second order perturbative contribution in J^n- In order to be theoretically 
consistent one should use perturbation theory up to the very same second 
order in J^n for all the quantities. E.g., pion spectrum is then determined 
by pion polarization operator expanded up to the very same order in f nN 

f-KN/ 

uj 2 ~ ml +fc 2 +77(f (w, k, q), I7*(lu, k, g) 

(21) 

The value 77° (u>, k, g) is easily calculated containing no any uncertain param- 
eters. For ui — > 0, k ~ pf of our interest and for isospin symmetric matter 

II* * -oo - ifto;, a c 2m »PF##» > o, /3b * > 0- (22) 

Replacing this value to (21) we obtain a solution with iw < already for 
£ > 0.3go that would mean appearance of the pion condensation. Indeed, 
the mean field begins to increase with the time passage ip ~ exp(— kot) ~ 
exp(at/p) until repulsive 7T7t interaction will not stop its growth. But it is 
experimentally proven that there is no pion condensation in atomic nuclei, i.e. 
even at g = g . The puzzle is solved as follows. FOPE model does not work 
for such densities. One should replace FOPE by the full TV TV interaction given 
by (4) , (7) . Essential part of this interaction is due to MOPE with vertices 
corrected by NN correlations. Also the TV TV -1 part of the pion polarization 
operator is corrected by NN correlations. Thus 

- i7 o fl (w,fc,0)7(<7V,fc,0) (23) 

being suppressed by the factor j(g',0J = 0, k ~ pp, g ~ g ) ~ 0.35 -j- 0.45. 
Final solution of the dispersion relation (21), now with full 77 instead of 77°, 
yields iw > for g = g whereas the solution with vjj < 0, which shows the 
beginning of pion condensation, appears only for g > g C7T > g . 



3 Neutrino cooling of neutron stars 
3.1 Emissivity of MMU process 

Since DU process is forbidden up to sufficiently high density g c u, the main 
contribution for g < g c u and T opac > T > T c comes from MMU pro- 
cesses schematically presented by two diagrams of Fig. 1. MNB reactions 
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give smaller contribution [23] . For densities g <§C go the main part of the NN 
interaction amplitude is given by the residual NN interaction. In this case 
the NN interaction amplitude can be better treated within the T matrix 
approach which sums up the ladder diagrams in the particle-particle chan- 
nel rather then by LM parameters. Calculations of MNB processes with the 
vacuum T matrix [77] found essentially smaller emissivity then that given 
by FOPE. Also the in-medium modifications of the T matrix additionally 
suppress the rates of both MMU and MNB processes, see [78] . Thus even at 
small densities the FOPE model may give only a rough estimate of the emis- 
sivity of two nucleon processes. At g > (0.5 -j- 0.7) go reactions in particle-hole 
channel and more specifically with participation of the soft pion mode begin 
to dominate. 

Evaluations [23,42,43,16] showed that the dominating contribution to 
MMU rate comes from the second diagram of Fig. 1, namely from contri- 



butions to it given by the first two diagrams of the series 




n n n n (24) 

whereas the third diagram, which naturally generalizes the corresponding 
MU(FOPE) contribution, gives only a small correction for g > g . The emis- 
sivity from the two first diagrams in a simplified notation [16,26] reads 



MMU 



[MOPE] 



au>„± (p Fr 



2.4 • 10 24 T 9 8 ( ± 
1 go 



r 8 FiC(A0 C(A> 



10/3 



«) 3 S 



to 



N 



au> v o(p Fn ) 



erg 



(25) 



where Tg = T/10 9 K is the temperature, to* and m* are the nonrelativistic 
effective neutron and proton masses, and the correlation factor r 8 is roughly 



r 8 - i}^ - PFe, q - PFeh 2 (g nn , u ~ o, k = p Fn )j 4 (g', o,p Fn ), 

2, x 7 2 (/>,9) + 3 5 i7 2 (5>,9) 
7 a{u,q) = 



and the second term in the factor 
F 1 ~ 1 



47 2 ( 5 ', 0,p F n)7j(^ - PFe, 9 - PFe) 



2/3 



(26) 



(27) 



is the contribution of the pion decay from intermediate states (first diagram 
(24)). The quantity r effectively accounts for an appropriate product of the 
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NN correlation factors in different irNiN 2 vertices. For charged pions the 
value fj, v ^ is incorporated in the expression for the effective pion gap, for 
neutral pions = 0. The value a ~ 1 depends on condensate structure, 
a = 1 for q < g C7T , and a = \[2 taking account of the new excitations on the 
ground of the charged 7r condensate vacuum for g > g C7r . The factor 

t(A ^ / cxp (-^/T) T<T cN , 

^ A ^-\l T > T c tv, N — (n, p) < 28 ) 

estimates the suppression caused by the nn and pp pairings. Deviation of 
these factors from simple exponents can be incorporated as in [35]. 

The ratio of the emissivities of MMU(MOPE) and MU(FOPE) is roughly 

i^lFOPE] " 10 ^ (p Fn )^ ± {p Fn ) ( ^ o) ' (29) 

For q ~ go this ratio is ~ 10 whereas being estimated with the only third 
diagram (24) it would be less then unit. 

3.2 Emissivity of DU-like processes 

NPBF processes. The one-nucleon processes with neutral currents given by 
the second diagram (18) for N = (n,p) are forbidden at T > T c by energy- 
momentum conservations but they can occur at T < T c . Then physically the 
processes relate to NPBF, see Fig. 2. However they need special techniques 
to be calculated [32,24]. These processes n — -> nvv and p — > pvv play very 
important role in the cooling of superfluid NS, see [24,25,16,26,36,35]. The 
emissivity for 3 types of neutrinos is given by [24,25] 

£[nPBF] = 3-4G 2 (^ 2 (fnn)+&gfj 2 (gnn))pFnm* n Al i ^A^ 



cm" • sec 



where £i = 1, £2 = for 5-pairing and £1 = 2/3, £2 = 4/3 for P-pairing, 
compare with the result [35] where £1 = 1, £2 = 2 were obtained. I removed 
some misprints existed in [32,24,25,16]. For neutrinos u>(q) = q and correla- 
tions are not so essential as it would be for uj <C q. Taking 7 2 ~ 1.3 in the 
range of S'-pairing we get C — 5 whereas for the P-pairing with g* A ~ 1.1 
we obtain C, ~ 9, in agreement with numerical evaluations [24] (£ ~ 7) used 
then within the cooling code in [26], I(x) = J Q °° ch 5 ydy/ (exp(xchy) + l) 2 , 
I(x 3> 1) — exp(— 2A/T)^y nT /AA, that serves an appropriate asymptotic 
(30). Emissivity of the process p — ► pvv is given by [24] 

£ (PPBF) = (±\ T< Tcp(31) 
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and £ 2 = for protons paired in S-state in NS matter, £3^1 is due to 
the second diagram (19) and has a complicated structure [73,76]. For the 
process (31) the part of NN and ee correlations is especially important. One 
has Kp p <~ 0.05 -7- 1 for 1 3qo and £3 <~ 1, and K pp + £3 ~ 1 instead of a 
small Cy ~ 0.006 value in absence of correlations, see discussion in subsection 
2.3. Thereby, in agreement with [24,26,73,76], the emissivity of the process 
p — > pvv can be compatible with that for n — ► n^i> in dependence on the 
relation between A p and A n . 

NPBF processes are very efficient for T < T c competing with MMU pro- 
cesses. The former win the content for not too massive stars. Analysis of 
above processes supports also our general conclusion on the crucial role of 
in-medium effects in the cooling scenario. 

Pion (kaon) condensate processes. The P wave pion condensate can 
be of three types: tt+ , 7t , and tt° with different values of the critical densities 
Qc-k = (g C 7r± 7 Q C7T + 1 QcTT a ), see [57]. Thus above the threshold density for the 
pion condensation of the given type, the neutrino emissivity of the MMU 
process (25) is to be supplemented by the corresponding PU processes 



X £ X 




The emissivity of the charged pion condensate processes with inclusion of the 
NN correlation effect (in a simplified treatment) renders, see [22,16], 



e[PU] ~ 7- 10 26 t£lL^P^ {pFe7PFe) ^(g',0,p Fn )T$ S m 2 6^^(33) 



Here g > g C7T and sin6> ~ ^/2\ui 2 \/m'^ for 6 <C 1. Of the same order of 
magnitude are emissivities of other possible ir condensate reactions, e.g. for 
7i7r" — > npev process at 9 <C 1 the numerical factor is about two times larger. 
Since tt^ condensation probably reduces the energy gaps of the superfluid 
states by an order of magnitude, see [79], we may assume that superfluidity 
vanishes above g C7T . Finally we note that though the PU processes have gen- 
uinely one-nucleon phase-space volumes, their contribution to the emissivity 
is suppressed relative to the DU by an additional j 2 (g' } 0,PFn) suppression 
factor due to existence of the extra (ttNN) vertex in the former case. 

Fig. 4 compares the mass dependence of the neutrino cooling rates L u /Cy, 
L v is the neutrino luminosity, Cy is the heat capacity, for MMU(MOPE) 
and MU(FOPE) for non-superfluid matter. For the solid curves, the neutrino 
emissivity in pion-condensed matter is taken into account according to (33) 
and (25) with the parameter a = \/~2. To be conservative we took q C7T — 
3qq. Dashed curves correspond to the model where no pion condensation is 
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allowed. As one sees, the medium polarization effects included in MMU may 
result in three order of magnitude increase of the cooling rate for the most 
massive stars. Even for stars of a rather low mass the cooling rate of MMU 
is still several times larger than for MU because even in this case the more 
efficient rate is given by the reactions shown by the right diagram in Fig. 1 
(first two diagrams of (24)). The cooling rates for the NS of M = 1.8M with 
and without pion condensate differ only moderately (by factor of 5 in this 
model). If we would take g C7T to be smaller the ratio of emissivity PU to MMU 
would decrease and could even become 1 in the vicinity of g c1T . Opposite, 
the reaction rates for the FOPE model are rather independent of the star's 
mass for the stars with masses below the critical value 1.63Mq, at which 
transition into the pion condensed phase occurs, and then jump to the typical 
PU value. It is to be stressed that contrary to FOPE model, the MOPE model 
[23] consistently takes into account the pion softening effects for g < g cv and 
both the pion condensation and pion softening effects on the ground of the 
condensate for g > p C7r . For g > g c x the kaon condensate processes come into 




0.5 1.0 1.5 2.0 

M/M, nB 

Fig. 4. Cooling rate due to neutrino emission as a function of star mass for a 
representative temperature of T = 3 x 10 s K. Superfluidity is neglected. 



play. Most popular is the idea of the S wave K~ condensation (e.g. see [13]) 
which is allowed at /i e > m* K _ due to possibility of the reaction e — > K~v. 
Analogous condition for the pions, [i e > to* _ , is not fulfilled owing to a strong 
S wave ttNN repulsion [57,16] (again in-mcdium effect!) otherwise S wave 
7r~ condensation would occur at smaller densities then K~ condensation. 
The neutrino emissivity of the K ~ condensate processes is given by equation 
analogous to (33) with a different NN correlation factor and an additional 
suppression factor due to a small contribution of the Cabibbo angle. However 
qualitatively scenario that permits kaon condensate processes is analogous to 
that with the pion condensate processes. 
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Other resonance processes. There are many other reaction channels 
allowed in the medium. E.g, any Fermi liquid permits propagation of zero 
sound excitations of different symmetry related to the pion and the quanta 
of a more local interaction determined here via f a ^ and g a ,p- These excita- 
tions being present at T 7^ may also participate in the neutrino reactions. 
The most essential contribution comes from the neutral current processes [23] 
given by first two diagrams of the series 

v 




Here the dotted line is zero sound quantum of appropriate symmetry. These 
are the resonance processes (second, of DU-type) analogous to those processes 
going on the condensates with the only difference that the rates of reactions 
with zero sounds are proportional to the thermal occupations of the corre- 
sponding spectrum branches whereas the rates of the condensate processes 
are proportional to the modulus squared of condensate mean field. Contribu- 
tion of the resonance reactions is as a role rather small due to a small phase 
space volume (q <~ T) associated with zero sounds. Please also bear in mind 
an analogy of the processes (34) with the corresponding phonon processes in 
the crust. 

DU processes. The proper DU processes in matter, as n — > pe~D e and 
pe~ — > nv e , 

e e 




v + I (35) 



P P 

should also be treated with the full vertices. They are forbidden up to the 
density g c u when triangle inequality pp n < pf p + PFe begins to fulfill. For 
traditional EQS like that given by the variational theory [18] DU processes 
are permitted for g > 5 go- The emissivity of the DU processes renders 

where m = fj, e = fi^ is the chemical potential of the leptons in MeV. In 
addition to the usually exploited result [15], (36) contains 7I pre-factor (15) 
due to NN correlations in the (3 decay vertices, see [24,26]. 

It was realized in [80] that the softening of the pion mode in dense neutron 
matter could also give rise to a rearrangement of single-fermion degrees of 
freedom due to violation of Pomeranchuk stability condition for g = g cF < 
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g C7r . It may result in organization of an extra Fermi sea for g c p < g < g C7r 
and at small momenta p < Q.2pp n , that in its turn opens a DU channel 
of neutrino cooling of NS from the corresponding layer. Due to the new 
feature of a temperature-dependent neutron effective mass, to* oc 1/T, we 
may anticipate extra essential enhancement of the DU process, corresponding 
to a reduction in the power of the temperature dependence from T 6 to T 5 . 
At early hot stage this layer becomes to be opaque for the neutrinos slowing 
the transport from the massive NS core to the exterior. 

Basing on the Brown-Rho scaling idea [71], we argued in [81] for the 
charged g meson condensation at a relevant density (g ce <~ 3go if m* would 
decrease to that density up to m s /2). If happened, it would open DU reaction 
already for g < g cg and close it for g > g cg due to an essential modification 
of the nuclear asymmetry energy. 

3.3 Comparison with soft X ray data 
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Fig. 5. Cooling of non-superfluid NS models of different masses constructed for 
the HV EOS [26]. Two graphs refer to cooling via MU(FOPE) +PU (left) and 
MMU(MOPE)+PU (right). In both cases, pion condensation is taken into account 
for the solid curves at g > q c -w, i.e. for M > 1.6Mq (q C7T is chosen to be 3^o) • The 
dashed curves in the right graph refer to a somewhat larger value of uj 2 (pFn), with- 
out pion condensation. The observed luminosities are labeled by dots. Possibilities 
of Fermi sea rearrangement and condensation are ignored. 

The heat transport within the crust establishes homogeneous density pro- 
file at times ^ (1 -f- 10)yr. After that time the subsequent cooling is deter- 
mined by simple relation CyT — —L, where Cy = J2i Cv,i and L = J2i Li 
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arc sums of the partial contributions to the heat capacity (specific heat in- 
tegrated over the volume) and the luminosity (emissivity integrated over the 
volume) . 

The nucleon pairing gaps are rather purely known. Therefore one may 
vary them. The "standard" and "nonstandard" scenarios of the cooling of NS 
of several selected masses for suppressed gaps are demonstrated in the left 
panel of Fig. 5, [26]. Depending on the star mass, the resulting photon lumi- 
nosities are basically either too high or too low compared to those given by 
observations. Situation changes, if the MMU process (25) is included. Now, 
the cooling rates vary smoothly with the star mass (see right panel of Fig. 
5) such that the gap between standard and non-standard cooling scenarios is 
washed out. More quantitatively, by means of varying the NS mass between 
(1 1.6) Mq, one achieves an agreement with a large number of observed 
data points. This is true for a wide range of choices of the lo 2 parameteriza- 
tion, independently whether pion (kaon) condensation can occur or not. Two 
parameterizations presented in Fig. 5 with pion condensation for g > 3qo 
and without it differ only in the range which is covered by the cooling curves. 
The only point which does not agree with the cooling curves belongs to the 
hottest pulsar PSR 1951+32. Other three points which according to Fig. 5, 
right panel, are also not fitted by the curves can be easily fitted by slight 
changes of the model parameters. The high luminosity of PSR 1951+32 may 
be due to internal heating processes, cf. [10]. 
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Fig. 6. Cooling of NS with different masses constructed for the HV EOS [26]. The 
cooling processes are MU-VS86, PU (only solid curves), NPBF, PPBF, and DU 
(only in the right graph). The dashed curves refer to the M — 1.7 and 1.9Mq 
models without pion condensate. Possibilities of Fermi sea rearrangement and g ± 
condensation are ignored. 



24 D.N. Voskresensky 



Wc turn now to cooling simulations where the MU, NPBF, DU and PU 
take place simultaneously. Parameters of the pairing gaps are from Fig. 6 of 
[26]. Fig. 6 shows the cooling tracks of stars of different masses, computed for 
the HV EQS. Very efficient at T < T c become to be NPBF processes which 
compete with MMU processes. The former prevail for not too massive stars in 
agreement with estimation [24] . The DU process is taken into account in the 
right graph, whereas it is neglected in the left graph. The solid curves refer 
again to the Cj 2 parameterization with phase transition to pion condensate, 
the dashed curves to the one without phase transition (see Fig. 4). For masses 
in the range between 1.0 and 1.6 M Q , the cooling curves pass through most of 
the data points. We again recognize a photon luminosity drop by more than 
two orders of magnitude for the 1.7 Mq mass star with pion condensate, due 
to suppression of the pairing gaps in this case. This drop is even larger if the 
DU is taken into account (right graph). This allows to account for the photon 
luminosity of PSR 1929+10. 

Thus, comparison with the observed luminosities shows that one gets 
quite good agreement between theory and observations if one includes into 
consideration all available in-medium effects assuming that the masses of the 
pulsars are different. We point out that the description of these effects is 
constructed in essentially the same manner for all the hadron systems as NS, 
atomic nuclei and heavy ion collisions, cf. [16,58]. 



3.4 Neutrino opacity 



Important part of in-medium effects for description of neutrino transport at 
initial stage of NS cooling was discussed in [23,42,43,16], where correlation ef- 
fects, pion softening and pion condensation, the latter for g > p C7r , were taken 
into account. The neutrino/antineutrino mean free paths can be evaluated 
from the corresponding kinetic equations via their widths F v (v) = —2lmII^-^ , 

where II R is the retarded self-energy, or within the QPA for the nucleons they 
can be also estimated via the squared matrix elements of the corresponding 
reactions. Thereby the processes which the most efficiently contribute to the 
emissivity are at early times (for T > 1 MeV) responsible for the opacity. 

In the above "nuclear medium cooling scenario" at T > T c the most es- 
sential contribution was from MMU. Taking into account of NN correlations 
in the strong coupling vertices of two-nucleon processes like MMU and MNB 
suppresses the rates, whereas the softening of the pion propagator essentially 
enhances them. For rather massive NS MOPE wins the competition. The 
mean free path of neutrino/antineurtino in MMU processes is determined 
from the same diagrams (24) as the emissivity. Its calculation (see (25)) with 
the two first diagrams yields 
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From the relation \ v ~ R follows evaluation of T opac . With the only first 
diagram we get a simple estimate 

Tr°^n^— (38 ) 
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For averaged value of the density g ~ Qq corresponding to a medium-heavy 
NS (< 1.4M Q ) with £o 2 {p Fn ) ~ 0.8m2 j ^ ~ 7 ~ (0.3 - 0.4) we get T opac ~ 
(1 -j- 1.5) MeV that is smaller then the value T opoc ~ 2 MeV estimated with 
FOPE [6]. For g ~ 2go that corresponds to a more massive NS we evaluate 
?opac — (0.3 -j- 0.5) MeV. Thus pion softening results in a substantial decrease 
of neutrino/antineutrino mean free paths and the value of T opac . 

The diffusion equation determines the characteristic time scale for the 
heat transport of neutrinos from the hot zone to the star surface to <~ 
R 2 Cycr~ 1 T~ 3 / X v (a is the Stefan-Boltzmann constant), whereas it follows 
that t ~10 min. for T ~ 10 MeV and g ~ g , uj 2 ~ 0.8™^, T ~ 0.4, 
m* N /rriN — 0.9, and to becomes as large as several hours for g ~ (2 4- 3)go- 
These estimates demonstrate that more massive NS cool down more slowly 
at T > T opac and faster at subsequent times then the less massive stars. 

Due to in- medium effects neutrino scattering cross section on the neutrons 
shown by the first diagram (18) requires the NN correlation factor 

l nv {^ q) = nr ^| . (39) 

as follows from (16). This results in a suppression of of the cross sections for 
u) < qvp n and in an enhancement for u) > qv Fn . Neutrino scattering cross 
sections on the protons are modified by 

2/ \ K 2 {fn P ,fnn,u,q) + ?>g 2 A il{g nn ,w,q) 
7 pl >,<7) = 1 + 3g 2 . ( 4 °) 

that results in the same order of magnitude correction as given by (39). 

Also there is a suppression of the vN scattering and MNB reaction rates 
for soft neutrinos (w ^ (3 -j- 6)T) due to multiple NN collisions 

i 
i 

(41) 



(LPM effect). One may estimate these effects simply multiplying squared 
matrix elements of the vN scattering and the MNB processes by the cor- 
responding suppression pre-factors [28]. Qualitatively one may use a general 
pre-factor Co(u)) = w 2 /[u 2 + T^], where i~jv is the nucleon width and the w is 
the energy of v or vv pair. In some works, see [82,83], correction factor, like 
Co, was suggested at an ansatz level. Actually one does not need any ansatze 
reductions. OTF, see [28], allows to calculate the rates using the exact sum 
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rule. The modification of the charged current processes due to LPM effect is 
unimportant since the corresponding value of uj is ~ pFe 3> -Tzv- 

The main physical result we discussed is that in medium the reaction rates 
are essentially modified. A suppression arises due to NN correlations (for 
lo <C kvFN) and infra-red pre-factors (coherence effects), and an enhancement 
due to the pion softening (and NN correlations for uj <~ q) and due to opening 
up of new efficient reaction channels. The pion softening demonstrates that 
already for densities g < go the nucleon system begins to feel that it may 
have 7r condensate phase transition for g > g c1T , although this g C7T value might 
be essentially larger than g or even not achieved. 



4 The rate of radiation from dense medium. OTF 

Perturbative diagramss are obviously irrelevant for calculation of in-medium 
processes and one should deal with dressed Green functions. The QPA for 
fermions is applicable if the fermion width is much less than all the typical 
energy scales essential in the problem (/> <C w c h). In calculation of the emis- 
sivities of vv reactions the minimal scale is uj c h — 6T, averaged vv energy for 
MNB reactions. For MMU uj c h — PFe- For radiation of soft quanta of fixed 
energy uj < T, u) c h — UJ- Within the QPA for fermions, the reaction rate with 
participation of the fermion and the boson is given by [22,23] 
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> 



(42) 



For equilibrium (T^O) system there is the exact relation 



(< (p\(p x >)(p) = iD ++ | <p c 



exp(f) 
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(43) 

where (< {p\<fii >)(p) means the Fourier component of the corresponding 
non-equilibrium Green function and (p c is the mean field. Thus the rate of 
the reaction is related to the boson spectral function Ab and the width (Is) 
being determined by the corresponding Dyson equation, see (9). Ab is the 
delta-function at the spectrum branches related to resonance processes, like 
zero sound. The poles associated with the upper branches do not contribute 
at small temperatures due to a tiny thermal population of those branches. 
There is also a contribution to ImD fi proportional to ImTT^ given by the 
particle-hole diagram. Within the QPA taking Im part means the cut of the 
diagram. Thus we show [23] that this contribution is the same as that could 
be calculated with the help of the squared matrix element of the two-nucleon 
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process 



4- 2 




(44) 



This is precisely what one could expect using optical theorem. Thus unlim- 
ited series of all possible diagrams with in-medium Green functions (see (18), 
(32), (34), (35)) together with two-fermion diagrams (as given by (24)) and 
multiplc-fermion diagrams (like (41)) would lead us to a double counting. The 
reason is that permitting the boson width effects (and beyond the QPA for 
fermions also permitting finite fermion widths) the difference between one- 
fcrmion, two-fermion and multiplc-fermion processes in medium is smeared 
out. All the states are allowed and participate in production and absorption 
processes. Staying with the QP picture for fermions, the easiest way to avoid 
mentioned double counting is to calculate the reaction rates according to 
(42), i.e. with the help of the diagrams of the DU-like type, which already 
include all the contributions of the two-nucleon origin. Multiple NN collision 
processes should be added separately. On the other hand, it is rather inconve- 
nient to explicitly treat all one-nucleon processes dealing with different spe- 
cific quanta instead of using of the full NN interaction amplitude whenever 
it is possible. Besides, as we have mentioned, consideration of open fermion 
legs is only possible within the QPA for fermions since Feynman technique 
is not applicable if Green functions of ingoing and outgoing fermions have 
widths (that in another language means possibility of additional processes) . 
Thus the idea came [24,28] to integrate over all in-medium states allowing all 
possible processes instead of specifying different special reaction channels. 

In [24,28,27] it was shown that OTF in terms of full non-equilibrium 
Green functions is an efficient tool to calculate the reaction rates including 
finite particle widths and other in-medium effects. Applying this approach, 
e.g., to the antineutrino-lepton (electron, fi~ meson, or neutrino) production 
[24] we can express the transition probability in a direct reaction in terms of 
the evolution operator S, 



where we presented explicitly the phase-space volume of vl states; lepton oc- 
cupations of given spin, n; , are put zero for v and v which are supposed to be 
radiated directly from the system (for T < T opac ). The bar denotes statistical 
averaging. The summation goes over complete set of all possible intermediate 
states {X} constrained by the energy-momentum conservation. It was also 
supposed that electrons/muons can be treated in the QPA, i.e. with zero 
widths. Then there is no need (although possible) to consider them in inter- 
mediate reaction states. Making use of the smallness of the weak coupling, 




(45) 
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we expand the evolution operator as S w 1 — i J^°° T {Vw{x) S nuc \(x)}dxo , 
where Vw is the Hamiltonian of the weak interaction in the interaction repre- 
sentation, S nuc i is the part of the S matrix corresponding to the nuclear inter- 
action, and T{...} is the chronological ordering operator. After substitution 
into (45) and averaging over the arbitrary non-equilibrium state of a nuclear 
system, there appear chronologically ordered (G~~), anti-chronologically or- 
dered (G ++ ) and disordered (G H and G h ) exact Green functions. 

In graphical form the general expression for the probability of the lepton 
(electron, muon, neutrino) and anti-neutrino production is as follows 

I v 



v I 
representing the sum of all closed diagrams ( — iiT h ) containing at least one 
( — V) exact Green function. The latter quantity is especially important. Var- 
ious contributions from {X} can be classified according to the number N of 
G h lines in the diagram 




dWlf _ (1 - ni )d 3 q e dqf 
dt 
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I v I 
This procedure suggested in [24] is actually very helpful especially if the 
QPA holds for the fermions. Then contributions of specific processes con- 
tained in a closed diagram can be made visible by cutting the diagrams 
over the (H — ), ( — h) lines. In the framework of the QPA for the fermions 
G ^ = 2ir\nFS(e + \i — e° — RcS R (e + fi,p)) {tlf are fermionic occupations, 
for equilibrium np = l/[exp((e — fj,p)/T) + 1]), and the cut eliminating the 
energy integral thus requires clear physical meaning. This way one establishes 
the correspondence between closed diagrams and usual Feynman amplitudes 
although in general case of finite fcrmion width the cut has only a sym- 
bolic meaning. Next advantage is that in the QPA any extra G h , since it 
is proportional to rip, brings a small (T/ef) 2 factor to the emissivity of the 
process. Dealing with small temperatures one can restrict by the diagrams of 
the lowest order in (G h G H ), not forbidden by energy-momentum conser- 
vations, putting T = in all G ++ and G Green functions. Each diagram 
in (46) represents a whole class of pcrturbative diagrams of any order in the 
interaction strength and in the number of loops. 

Proceeding further we may explicitly decompose the first term in (46) as 
v I v I 




•XX XX 



(47) 



I 

The full vertex in the diagram (47) of given sign is irreducible with respect to 
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the (H — ) and (— +) nucleon-nucleon hole lines. This means it contains only 

the lines of given sign, all ( ) or (H — h). Second diagram with anomalous 

Green functions exists only for systems with pairing. In the framework of 
the QPA namely these diagrams determine the proper DU and also NPBF 
processes calculated in [24,25] within OTF. In the QP picture the contribution 
to DU process vanishes for g < g c u . Then the second term (46) comes into 
play which within the same QP picture contains two-nucleon processes with 
one (G h G H ) loop in intermediate states, etc. 

The full set of diagrams for 77 h can be further explicitly decomposed as 
series [28] (from now in brief notations) 

• GO • x - : -'OSO- * ... 

Full dot is the weak coupling vertex including all the diagrams of one sign, 
NN interaction block is the full block also of one sign diagrams. The lines 
are full Green functions with the widths. The most essential term is the one- 
loop diagram (see (47)), which is positive definite, and including the fermion 
width corresponds to the first term of the classical Langevin result, for details 
see [28]. Other diagrams represent interference terms due to rescatterings. In 
some simplified representations (e.g., as we used within Fermi liquid theory) 
the 4-point functions (blocks of NN interaction of given sign diagrams) be- 
have like intermediate bosons (e.g. zero-sounds and dressed pions). In general 
it is not necessary to consider different quanta dealing instead with the full 
NN interaction (all diagrams of given sign). For particle propagation in an 
external field, e.g. infinitely heavy scattering centers (proper LPM effect), 
only the one-loop diagram remains, since one deals then with a genuine one- 
body problem. In the quasiclassical limit for fermions (with small occupations 
rip) all the diagrams given by first line of series (48) with arbitrary number 

of " h" lines are summed up leading to the diffusion result, for details 

see [28] . For small momenta q this leads to a suppression factor of the form 
C = uj 2 / (uj 2 + T 2 ) , r x incorporates rescattering processes. In general case 
the total radiation rate is obtained by summation of all diagrams in (48). 
The value — i77 h determines the gain term in the generalized kinetic equa- 
tion for G h , see [84,85], that allows to use this method in non-equilibrium 
problems, like for description of neutrino transport in semi-transparent region 
of the neutrino-sphere of supernovae, as we may expect. 

In the QP limit diagrams 1, 2, 4 and 5 of (48) correspond to the MMU 
and MNB processes related to a single in-medium scattering of two fermionic 
QP and can be symbolically expressed as Feynman amplitude (49a) 
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(a) 



(c) 



I 





(49) 



The one-loop diagram in (48) is particular, since its QP approximant in many 
cases vanishes as we have mentioned. However the full one-loop includes 
QP graphs of the type (49b), which survive to the same order in r^/^ch 
as the other diagrams (therefore in [24] where QP picture was used this 
diagram was considered as allowed diagram). In QP series such a term is 
included in second diagram of (46) although beyond the QPA it is included 
as the proper self-energy insertion to the one-loop result, i.e. in first term 
(46) [28] . In fact it is positive definite and corresponds to the absolute square 
of the amplitude (49a). The other diagrams 2, 4 and 5 of (48) describe the 
interference of amplitude (49a) cither with those amplitudes where the weak 
coupling quantum (W pair) couples to another leg or with one of the exchange 
diagrams. For neutral interactions diagram (48:2) is more important than 
diagram 4 while this behavior reverses for charge exchange interactions (the 
latter is important, e.g., for gluon radiation from quarks in QCD transport 
due to color exchange interactions) . Diagrams like 3 describe the interference 
terms due to further rescatterings of the source fermion with others as shown 
by (49c). Diagram (48:6) describes the production from intermediate states 
and relates to the Feynman graph (49d). For photons in the soft limit (ui <C 
Ef) this diagram (49d) gives a smaller contribution to the photon production 
rate than the diagram (49a), where the normal bremsstrahlung contribution 
diverges like l/tu compared to the 1/e^-value typical for the coupling to 
intermediate fermion lines. For vv bremsstrahlung (49d) gives zero due to 
symmetry. However in some cases the process (49d) might be very important 
even in the soft limit. This is indeed the case for the MMU process considered 
above. Some of the diagrams which are not presented explicitly in (48) give 
more than two pieces, if being cut, so they never reduce to the Feynman 
amplitudes. However in the QPA they give zero contribution [28]. 

With Fp ~ ir 2 F 2 /ep for Fermi liquids, the criterion rp <C uj c h ~ T is 
satisfied for all thermal excitations As ~ T <C ep/ir 2 . However with the 
application to soft radiation this concept is no longer justified. Indeed series 
of QP diagrams is not convergent in the soft limit and there is no hope 
to ever recover a reliable result by a finite number of QP diagrams for the 
production of soft quanta. With full Green functions, however, one obtains 
a description that uniformly covers both the soft (w <C Tp) and the hard 
(uj 3> rp) regimes. In the vicinity of g C7T the quantity TV being roughly 
estimated in [62,64] as rp cx 7T 2 r 2 Tm 7T /oj, and coherence effects come into 
play. 
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In order to correct QP evaluations of different diagrams by the fermion 
width effects for soft radiating quanta one can simply multiply the QP results 
by different pre-factors [28]. E.g., comparing the one-loop result at non-zero 

with the first non-zero diagram in the QPA (Tp = in the fermion Green 
functions) we get 
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at small momentum q. For the next order diagrams we have 
+ „ _+ 

= Ci(w) 



-oo- 



C (w) 




QPA 



QPA 



2\2 ' 



(51) 



where factors C , Ci, ... cure the defect of the QPA for soft u). The factor C 
complies with the replacement u> — > lu + i-Tp . A similar factor is observed in 
the diffusion result, where however the macroscopic relaxation rate F x enters, 
due to the resummation of all rescattering processes. 

Finally, we demonstrated how to calculate the rates of different reactions 
in dense equilibrium and non-equilibrium matter and compared the results 
derived in closed diagram technique with those obtained in the standard 
technique of computing of the squared matrix elements. 



5 Cooling of Hybrid Neutron Stars 

Let us in addition review how the color superconducting quark matter, if 
exists in interiors of massive NS, may affect the neutrino cooling of HNS, 
[53] . A detailed discussion of the neutrino emissivity of quark matter without 
taking into account of the possibility of the color superconductivity has been 
given first in ref. [86]. In this work the quark direct Urea (QDU) reactions 
d — ► uev and ue — ► dv have been suggested as the most efficient processes. 
In the color superconducting matter the corresponding expression for the 
emissivity modifies as 

e Q° u ~ 9.4 x 10 26 a s (g/g a )Y c 1 / 3 C QDV T 9 6 erg cm" 3 s" 1 , (52) 

where due to the pairing the emissivity of QDU processes is suppressed 
by a factor, very roughly given by Cqdu ~ exp(— A q /T). At g/g — 2 
the strong coupling constant is a s w 1 decreasing logarithmically at still 
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higher densities, Y e = g e /g is the electron fraction. If for somewhat larger 
density the electron fraction was too small (Y e < Y ec ~ 10~ 8 ), then all 
the QDU processes would be completely switched off [87] and the neutrino 
emission would be governed by two-quark reactions like the quark modified 
Urea (QMU) and the quark bremsstrahlung (QB) processes dq — ► uqev and 
<7i<72 — * qxqivv, respectively. The emissivities of QMU and QB processes are 
smaller than that for QDU being suppressed by factor Cqmu ~ exp(— 2A q /T) 
for T < T cr it,<j — 0.4 A q . For T > T cr i t . q all the £ factors are equal to unity. 
The modification of T CT n iq (A q ) relative to the standard BCS formula is due 
to the formation of correlations as, e.g., instanton- anti-instanton molecules. 
The contribution of the reaction ee — ► eevv to the cmissivity is very small 
[88], e e v e ~ 10 12 {q/ QqY^T^ erg cm -3 s _1 , but it can become important 
when quark processes are blocked out for large values of A q /T in supercon- 
ducting quark matter. 

For the quark specific heat [53] used expression of [86] being however 
suppressed by the corresponding ( factor due to color superfluidity. Therefore 
gluon-photon and electron contributions play important role. 

The heat conductivity of the matter is the sum of partial contributions 
k = J2i K i> K T 1 — J2j K ij 1 i where i,j denote the components (particle 
species). For quark matter k is the sum of the partial conductivities of the 
electron, quark and gluon components k = n e + n q + n g , where n e ~ n ee is 
determined by electron-electron scattering processes since in superconducting 
quark matter the partial contribution 1/ ' K eq (as well as 1/ n gq ) is additionally 
suppressed by a Cqdu factor, as for the scattering on impurities in metallic 
superconductors. Due to very small resulting value of k the typical time 
necessary for the heat to reach the star surface is large, delaying the cooling 
of HNS. 

EQS used in ref. [53] included a model of the hadron matter, region of 
mixed phase and of pure quark matter. A hard EQS for the hadron matter 
was used, finite size effects were disregarded in description of the mixed phase, 
and the bag constant B was taken to be rather small that led to the presence 
of a wide region of mixed and quark phases already for the HNS of the mass 
M = 1.4 Mq. On the other hand, absence of a dense hadron region within 
this EQS allowed to diminish uncertainties in description of in-medium effects 
in hadronic matter suppressing them according to that used in the "standard 
scenario" . 

With the above inputs ref. [53] solved the evolution equation for the tem- 
perature profile. In order to demonstrate the influence of the size of the 
diquark and nucleon pairing gaps on the evolution of the temperature pro- 
file solutions were performed with different values of the quark and nucleon 
gaps. Comparison of the cooling evolution (lg T s vs. lg t) of HNS of the 
mass M = 1.4 M Q is given in Fig. 7. The curves for A q > I MeV are very 
close to each other demonstrating typical large gap behaviour. The behaviour 
of the cooling curve for t < 50 -j- 100 yr is in straight correspondence with 
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the heat transport processes. The subsequent time evolution is governed by 
the processes in the hadronic shell and by a delayed transport within the 
quark core with a dramatically suppressed neutrino emissivity from the color 
superconducting region. In order to demonstrate this feature a calculation 
was performed with the nucleon gaps (Ai(n), i = n,p) being artificially 
suppressed by a factor 0.2. Then up to lg(t[yr]) 4 the behaviour of the 
cooling curve is analogous to the one would be obtained for pure hadronic 
matter. The curves labelled "MMU" show the cooling of hadron matter with 
inclusion of appropriate medium modifications in the NN interaction. These 




iog(t[yr]) 

Fig. 7. Evolution of the surface temperature T s of HNS with M = 1.4M for 
T a = (10T m ) 2/3 , where T is in K, see [8]. Data are from [26] (full symbols) and 
from [35] (empty symbols), t w is the typical time which is necessary for the cooling 
wave to pass through the crust. 

effects have an influence on the cooling evolution only for lg(f [yr]) ^5 2 since 
the specific model EQS used does not allow for high nucleon densities in the 
hadron phase at given example of HNS of M = 1.4 Mq. The effect would 
be much more pronounced for larger star masses, a softer EQS for hadron 
matter and smaller values of the gaps in the hadronic phase. Besides, incor- 
poration of finite size effects within description of the mixed phase reducing 
its region should enlarge the size of the pure hadron phase. 

The unique asymptotic behaviour at lg(f[yr]) > 5 for all the curves cor- 
responding to finite values of the quark and nucleon gaps is due to a com- 
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petition between normal electron contribution to the specific heat and the 
photon cmissivity from the surface since small exponentials switch off all the 
processes related to paired particles. This tail is very sensitive to the inter- 
polation law T s — f(T m ) used to simplify the consideration of the crust. The 

2/3 

curves coincide at large times due to the uniquely chosen relation T s oc T„{ . 

The curves for A q =0.1 MeV demonstrate an intermediate cooling be- 
haviour between those for A q — 50 MeV and A q = 0. Heat transport becomes 
not efficient after first 5 4- 10 yr. The subsequent 10 4 yr evolution is governed 
by QDU processes and quark specific heat being only moderately suppressed 
by the gaps and by the rates of NPBF processes in the hadronic matter (up 
to lg(i[yr]) < 2.5). At lg(i[yr]) > 4 begins the photon cooling era. 

The curves for normal quark matter (A q = 0) are governed by the heat 
transport at times t 5 yr and then by QDU processes and the quark specific 
heat. The NPBF processes are important up to lg(£[yr]) < 2, the photon era 
is delayed up to lg(i[yr]) > 7. For times smaller than t w (see Fig. 7) the heat 
transport is delayed within the crust area [15]. Since for simplicity this delay 
was disregarded in the heat transport treatment, for such small times the 
curves should be interpreted as the T m {t) dependence scaled to guide the eye 
by the same law oc T r ?/ 3 , as T s . 

For the CFL phase with large quark gap, which expected to exhibit the 
most prominent manifestations of color superconductivity in HNS and QCNS, 
[53] thus demonstrated an essential delay of the cooling during the first 50 -j- 
300 yr (the latter for QCNS) due to a dramatic suppression of the heat 
conductivity in the quark matter region. This delay makes the cooling of 
HNS and QCNS not as rapid as one could expect when ignoring the heat 
transport. In HNS compared to QCNS (large gaps) there is an additional 
delay of the subsequent cooling evolution which comes from the processes in 
pure hadronic matter. 

In spite of that we find still too fast cooling for those objects compared 
to ordinary NS. Therefore, with the CFL phase of large quark gap it seems 
rather difficult to explain the majority of the presently known data both in 
the cases of the HNS and QCNS, whereas in the case of pure hadronic stars 
the available data are much better fitted even within the same simplified 
model for the hadronic matter. For 2SC (3SC) phases one may expect anal- 
ogous behaviour to that demonstrated by A q = since QDU processes on 
unpaired quarks are then allowed, resulting in a fast cooling. It is however 
not excluded that new observations may lead to lower surface temperatures 
for some supernova remnants and will be better consistent with the model 
which also needs further improvements. On the other hand, if future observa- 
tions will show very large temperatures for young compact stars they could 
be interpreted as a manifestation of large gap color superconductivity in the 
interiors of these objects. 

Concluding, the "nuclear medium cooling scenario" allows easily to 
achieve agreement with existing data. However there remains essential uncer- 



Medium Effects 35 



tainty in quantitative predictions due to a poor knowledge, especially, of the 
residual interaction treated above in an economical way within a phcnomeno- 
logical Fermi liquid model which needs further essential improvements. As 
for color superconductivity in HNS, QCNS, and QS, characterizing by large 
diquark pairing gaps, we did not find an appropriate fit of existing X ray 
data. Situation will be changed if more cold, old objects will be observed. 
Also, if young hot objects ( on scale ~ 10 2 yr.) will be observed it could be 
interpreted as a signature of CFL phase in these objects. 
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